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Abstract 

The quantum mechanics of position measurement of a macroscopic object is typically inaccessible 
due to strong environmental coupling and classical noise. Here we show that a micromechanical mem- 
brane resonator subject to an increasingly strong continuous position measurement exhibits a quantum 
mechanical back action force that rises in accordance with the Heisenberg uncertainty relation. For our 
optically-based position measurements, the specific form of the back action is that of a fluctuating radia- 
tion pressure from the Poisson-distributed photons in the coherent measurement field, termed radiation 
pressure shot noise. We demonstrate a radiation pressure shot noise drive that is comparable in magni- 
tude to the thermal forces in our system. Additionally, we observe a temporal correlation between the 
fluctuations in radiation force and position of the resonator that may also be interpreted as a quantum 
non-demolition measurement of the photon field. 

It is a well-known principle of quantum mechanics that any position measurement must be accompanied 
by a back action force that perturbs the momentum of the system. For an optically-based ideal position 
measurement, the manifestation of this back action force is termed radiation pressure shot noise (RPSN) 
[IJ|2j[3]. Here, for example, the Poisson distributed photons (shot noise) of a coherent optical field reflect 
off the object to be measured, transferring recoil momentum at random times. For a harmonically bound 
object measured over many oscillation periods and by the accumulation of many photons, the increased 
spread in momentum will lead to a corresponding increase in displacement fluctuations, obscuring the initial 
position. Additionally, shot noise directly sets an imprecision noise floor on any method of optical detection. 
The point in an ideal measurement where the imprecision from the sum of these two effects is minimized is 
termed the standard quantum limit (SQL) [3, 4|[5j|6]. 

This idealized picture of quantum measurement and back action for large objects is typically obscured 
by classical noise sources such as thermal motion. Indeed, no solid object has yet been put into the regime 
where RPSN is the dominant driving force. The back action effect of RPSN is more readily observed if 
the measurement strength is increased beyond the SQL. As illustrated in Fig. 1A, the level of back action 
increases linearly with optical power or measurement strength. Radiation pressure effects are expected to 
soon limit precision interferometric measurements such as the LIGO gravitational wave detector. In the next 
generation of LIGO, extremely high laser powers will enhance strain sensitivity at high frequencies; with 
such powers RPSN is predicted to dominate the sensitivity at intermediate frequencies, even with large test 
masses at room temperature [2j [7] . Such quantum limits to continuous position measurement have long been 
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foreseen, and many ideas have been developed to circumvent the SQL by employing quadrature-squeezed 
light [8] or back action evasion techniques [9 . 

While a shot noise imprecision noise floor is a ubiquitous measurement limitation, the experimental 
signatures of RPSN have remained elusive. Effects of RPSN have been observed in samples of dilute, laser 
cooled atomic gas [TO] [11] , and several groups have made both experimental and theoretical studies of the 
predicted signatures of RPSN [12] [4j [I3j [14] [15] [16] . A new promising direction is low-mass, high-frequency, 
and high-optomechanical-coupling mechanical resonators that have been developed in the context of cavity 
optomechanics. Their enhanced response to optical forces makes them an ideal platform to study the effects of 
RPSN. In such devices analogous effects of electron shot noise measurement back action have been observed 
[T7] . Also, measurements of a weakly optically probed nanomechanical resonator have been interpreted 
as arising from correlations between sensing noise and back action [18]. Some of these nanomechanical 
resonators have even recently been cooled with electromagnetic radiation to near their quantum mechanical 
ground state p~9j [20] [21] , illustrating the capacity for dominant coherent optical forces. In this report, we 
demonstrate a micro-optomechanical resonator where the effects of RPSN are comparable in magnitude to 
those of the thermal forces. We observe the increased average amplitude of motion expected from RPSN; we 
also demonstrate temporal correlations between the optical driving force and the resonator position that can 
be interpreted as a quantum non-demolition (QND) measurement of intracavity photon fluctuations using 
the mechanical object as a meter. 

Our optomechanical system consists of a silicon nitride membrane microresonator [22] [23] inside of a 
Fabry-Perot optical cavity specially designed to operate at cryogenic temperatures (See Fig. [2] and Ref. 
[24]). Pioneering work by a group at Yale [22], has shown that such systems can achieve both large op- 
tomechanical coupling and very weak environmental coupling through optical and mechanical loss channels. 
Optomechanical coupling is achieved through a dispersive interaction, where the cavity resonance frequency 
shifts with the location of membrane along the optical standing wave. This interaction imprints phase and 
amplitude modulation on laser light transmitted through the system, which allows for measurement of the 
membrane motion. In conjunction, the laser applies an optical gradient force to the membrane, pushing 
it toward higher optical intensity. The mechanical modes are those of a highly tensioned square plate 0.5 
mm on a side and 40 nm thick, with an effective mass of about 7 ng. We operate in a helium flow cryostat 
at a temperature of 4.9 K, where intrinsic mechanical linewidths, ro/27r, are typically less than 1 Hz. For 
the (2,2) mode oscillating at cj m /27r = 1.55 MHz (twice the fundamental mode frequency), we achieve a 
maximum single photon coupling rate g/2it = 17 Hz. 

We use two laser fields in the cavity, derived from the same 1064 nm laser source, both coupled to the 
same spatial mode of the cavity, but with orthogonal polarizations [I2j [13] . All other optical modes are far 
enough separated in frequency to be neglected. The half-planar, 5.1 mm length, 72 fim mode waist cavity 
typically has a full linewidth k/2tt ~ 1 MHz, which varies slightly with the membrane position. The "signal" 
beam, typically of much larger intensity than the other, is maintained on the optical resonance via active 
feedback. This beam provides the RPSN. Its transmitted shot noise intensity fluctuations constitute a record 
of the optical force on the resonator. The corresponding sensitive position measurement is wholly imprinted 
in the unrecorded phase modulation spectrum. Additional phase noise from fluctuations in the cavity-laser 
detuning would preclude fully shot-noise-limited phase-quadrature detection in this system [24]. The much 
weaker "meter" beam is tuned to the red of the optical resonance. The resonator's displacement spectrum is 
imprinted in the transmitted intensity spectrum of this laser. While its shot noise drive is much smaller, the 
meter beam provides optical Raman sideband cooling of the mechanical resonator [25] . The optical damping 
greatly eases the requirements on the relative detuning of the signal beam from the cavity resonance due 
to both the onset of parametric instabilities at positive detuning and the contamination of cross correlation 
measurements by thermal motion near the mechanical resonance [14] [16]. However, Raman cooling does not 
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Figure 1: Displacement spectrum measurements. (A) Canonical picture of continuous position measurement. 
Contributions are: RPSN (black), thermal motion (orange), and zero point motion (purple). They combine 
to give the expected result of our measurement (blue). The dashed curve represents the effective displacement 
noise from the shot noise limited imprecision of an ideal measurement. (B) Spectral density for signal beam 
photocurrent near uj m (blue circles) for device A and measured shot noise floor (black). (C) Transmission 
spectra for device B at Ns = 2.0 x 10 T (blue) and Ns = 3.6 x 10 8 (orange). (D) Plotted versus signal 
power P are: Peak displacement spectral density for device A (green diamonds) and device B (red squares), 
expected values (blue), thermal contribution (orange), bounds on theoretical estimates including systematic 
uncertainty in device parameters and classical noise contribution (dashed), and expected RPSN contribution 
(black). Parameters for Device A (Device B): g/2ir = 16.3 ± 0.6 Hz (16.1 ± 0.3 Hz), k/2tt = 1.17 MHz 
(0.89 MHz), A s /2tt = 1.5 ± 1.5 kHz (2 ± 0.5 kHz), A m /2tt = 1.6 MHz (0.7 MHz), N M = 3.4 ± 0.3 x 10 6 
(7.0 ± 0.3 x 10 6 ), A^ ax = 1.2 ± 0.1 x 10 8 (4.4 ± 0.1 x 10 8 ), cj m /27r = 1.575 MHz (1.551 MHz), r /27r = 0.116 
Hz (0.47 Hz), r m /27r = 3 kHz (1.43 kHz). 
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increase the ratio of peak displacement due to RPSN versus the thermal displacement, as the resonator's 
response to both forces is damped equally in the cooling process. 




■►I 

Signal 



Meter 



Mirror Mirror 
4 He Cryostat 



Figure 2: Experimental Setup. The signal and meter beams are combined and separated with polarizing 
beam splitters (PBS) before and after interacting with the optomechanical cavity, and detected directly on 
photodetectors (PD). The inset photograph shows an in-situ image of the square membrane and optical 
mode spot, with dashed lines indicating the nodes of (2,2) mechanical mode. The plot on the right shows 
spectra from device B of the photocurrents Si s (uj)/Ig (green), Si m (uj)/I m (red), as well as the noise floors 
that include detector noise and the dominant shot noise (gray). 



The effect of the optomechanical coupling on the resonator from a single laser [25 , 26 , or multiple beams in 
orthogonal modes [14] has been well studied. The resonator's mechanical susceptibility is modified to include 
optomechanical damping and frequency shifts from each laser. Additionally, the amplitude of motion, or 
effective phonon occupation is modified. The optomechanical damping cools the resonator, while RPSN 
increases the amplitude of motion. In equilibrium, an effective phonon occupation can be computed from a 
simple rate equation as n m = (n t h^o + TisTs + tlm^ m) /^m- Here n t h is the thermal phonon occupation; 
ns and Ts (jim and Tm) are the effective bath temperature and optomechanical damping rate of the signal 
(meter) laser. The total mechanical damping rate is T m = r + Ts In our experiments Tm ^> r , T^, 

while As ~ and Ns ^> Nm, where A 5, Ns (Am, Nm) are the detuning from cavity resonance and 
intracavity photon occupation of the signal (meter) beam. The RPSN dominates over thermal noise when 
Rs — Cs/ n t\iO- + (2c<J m /ft) 2 ) -1 > 1, where Cs = 4:Nsg 2 / kTq is the multiphoton cooperativity. We are 
able to reach this high cooperativity regime due to the small mass, weak intrinsic damping, and cryogenic 
environment of our resonator. 

The increase in phonon occupation due to RPSN is shown in Fig. [I] Transmission spectra, S/ M (u;), for 
the meter beam are shown in Fig. [Tj3 for the case of negligible Ns and for Ns sufficiently large that Rs ~ 1. 
We observe a marked increase in the spectral area and hence n m . Fig. [lp shows, for two different devices, 
A and B, the peak displacement spectral density, S z (oj m ), inferred from transmitted intensity spectra. The 
parameters of each device are similar, except To is much smaller for device A, while device B has a much lower 
optical loss. Thus, device A enters the RPSN dominated regime at lower Ns than device B, while device B 
provides higher signal-to-noise measurements. We scale the peak displacement spectral density by the SQL 
peak spectral density S^ L (uj m ) = fi/muJmTm [5], at the damping provided by the meter beam used in the 
measurement. The optical measurement strength is scaled such that p/p s Q L = 1 at the SQL. Working in 
these quantum measurement units allows for direct comparison of the two devices. The increased spectral 
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density also includes a small contribution from classical radiation pressure noise. The total transmitted noise 
spectral density of the signal beam is shown in Fig[lp. For device A, at all measured laser powers, the noise 
level is within a few percent of the measured shot noise level at the photodetector. This noise level, referred 
back to the intracavity level, puts an upper bound of at most 15% of the displacement signal arising from 
classical laser noise. Taking into account the thermal motion and classical intensity noise, we can attribute at 
least 55% of the total displacement spectrum to RPSN for device A at the maximum signal beam strength. 
For device B this lower bound is 40%. The dashed curves of of Fig. ID represent bounds on the expected 
spectral densities accounting for systematic uncertainties in the device parameters including uncertainty in 
the classical noise level. Another effect that might mimic the effect of RPSN is physical heating. To test for 
physical heating, we monitor the response of a higher frequency, weakly optomechanically coupled mechanical 
mode where RSPN effects are expected to be negligible. This response indicates the bath temperature for 
the higher-frequency mode, and hence for our mode of interest if we assume a common thermal bath, does 
not rise by more than 10% due to absorbed laser light. Details of the calibration process are given in the 
Supplemental Material. 

We now turn to the study of the temporal correlations between the signal beam photocurrent, which 
records the radiation pressure drive, and the meter beam photocurrent, which contains the mechanical 
response to the RPSN. We compute the spectrum of the two-time cross correlation function Si sm (uj) = 
Im(w)), where I(uu) is the complex Fourier transformation of the photocurrent I(t), and the angle 
brackets represent an average over many realizations of the experiment. In the limit T m <C ft, the correlation 
should reflect the Lorentzian response function of the optically damped harmonic oscillator, driven by an 
approximately white noise source. Fig. [3]A. shows a cross correlation measurement for device B, and for 
reference, Si s (uj) x Si m (uj) is shown in the same plot, and the individual traces of Si s and *S/ M for these 
data are shown in Fig. |2j For a resonator well into the RPSN dominated regime, measured with ideal 
quantum efficiency, the two curves should coincide at the peak, yielding a perfect correlation. However, a 
white, uncorrelated measurement noise background, dominated by the meter beam's shot noise, separates 
the two curves away from the mechanical resonance. Additionally, the peak value of the cross correlation 
is smaller, arising mainly from two sources. The residual thermal motion is uncorrelated with the radiation 
pressure drive. Also, the imperfect detection efficiency of the signal laser leads to a loss of correlation, as 
information about the intracavity shot noise drive is unrecorded. We measure a peak relative correlation 
C(u m ) = \S IsM (uj m )\ 2 /S Is (uJm)Si M (uJm) = 0.14. A simple estimate is C(u m ) = Rs/0- + Rs) x k r /k x ei x 
^darki = 0.15 ± 0.02, where Rs/(1 + Rs) — 0.4 ± 0.03 is the fraction of the displacement spectrum due to 
RPSN, kr/k — 0.59 is the fraction of the light exiting the cavity through the output port, e\ — 0.79 ± 0.03 
is the post-cavity detection efficiency, and edarki = 0.8 is due to dark noise on the photodetector. By 
intentionally adding white classical intensity noise of much larger amplitude than shot noise to the signal 
laser, we are able to demonstrate in Fig. [3^3 a (classical) relative cross correlation that approaches unity. 

A detailed analysis of the cross correlation response yields considerable information about the system. 
Fig. |3p shows the phase of the correlation both with and without large classical intensity noise on the 
signal beam. Both show the 180 degree phase shift expected from the mechanical response. However, we 
also expect a phase offset of arctan (2uo rn / n) between the classical noise dominated drive and the shot noise 
dominated drive (See Supplemental Material and Ref. [H]). Fitting these two curves yields a relative phase 
shift that is slightly smaller than predicted, implying that in our measurement with no intentionally added 
classical noise, about 75 % of the radiation pressure drive is from shot noise, in agreement with estimates 
obtained by directly measuring the classical noise level of the signal beam transmission spectrum. 

If the signal beam is not precisely on resonance then the cross correlation will be distorted. Mechanical 
motion transduced directly onto the signal beam may either constructively or destructively add to the 
RPSN correlation depending on the sign of A5. By fitting the correlation data to the expected lineshape (see 
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Supplemental Material) with As as a fit parameter, we estimate As — 0.0003ft, which implies a contribution 
of only 3 % to the peak cross correlation from thermal motion. Also, note that since the thermal contribution 
is linear in As, the effect of the few kilohertz residual technical fluctuations in the laser-cavity detuning should 
average to that of the mean value. We have also performed an experimental test to demonstrate the rejection 
of ambient motion from the cross correlation spectrum. Here, we mechanically excite the membrane with 
a white- noise driven piezoelectric actuator. As shown in Fig. [3^3, the membrane is excited to a greater 
amplitude of motion near the mechanical resonance frequency. However, the cross correlation spectrum 
remains unchanged, implying very little of the ambient motion is transduced. 
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Figure 3: Cross correlation measurements. (A) The magnitude squared of Si SM {uu) / IsIm for device B (red), 
expected cross correlation including systematic uncertainty (gray), Si s (uo) x Si m (lj) (black). Parameters are 
as listed in Fig. [l] caption except A s /2ir = 300 ± 100 Hz, g/2ir = 14.8 ±0.4 Hz, N s = 3.2 x 10 8 (B) S Ism (uj) 
(green) and, Si s (u) x Si M (w) (orange) where classical intensity noise at the level of ~ 40 times shot noise, is 
added to the signal beam, Si SM (u) (blue) Si s (u) x Si M (uo) (purple) where membrane is driven with excess 
mechanical noise. Fits to the data of part (A) are displayed for reference (dashed black). (C) Phase of the 
cross correlation with classical intensity noise on signal beam (green) and without (red). Black curves are 
fits to the data. 



An alternative interpretation of these cross correlation measurements is that of a QND measurement of 
the intracavity photon fluctuations of the signal beam [27]. Here, the membrane acts as the measurement 
device, with its state of motion recording the photon fluctuations over the band of the mechanical resonance. 
In fact our C parameter is equivalent to the state preparation fidelity developed in Ref. [28] for nonideal 
QND measurements. Further, it has been shown that frequency-dependent ponderomotive squeezing of the 
signal beam quantum noise is possible [29] , and such squeezing has recently been demonstrated in an atomic 
gas cavity optomechanical system [30]. However for the laser configuration used in these measurements 
(As = 0), we do not expect to see squeezing of the detected amplitude quadrature. In conclusion, we have 
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made several measurements of the effects of radiation pressure shot noise on a micromechanical resonator. 
These effects are signatures of quantum measurement back action from continuous position measurement 
by an optical probe. Our observations open the door to realizing position measurement near the SQL if 
residual thermal noise and excess cavity-laser phase noise can be eliminated. We have also demonstrated 
the parameters necessary for strong ponderomotive optical squeezing. 
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Supplementary Materials for: 
"Observation of Radiation Pressure Shot Noise" 

T. P. Purdy, R. W. Peterson, and C. A. Regal 

S.l Theory 

Here we calculate the response of an optomechanical system to two independent laser driving fields. Addi- 
tionally, we compute the spectrum of the two time cross correlation of the photocurrents of the transmitted 
laser fields. This formalism is used to model our membrane cavity optomechanical system, as well as our 
specific detection setup that incorporates direct photodetection of the transmitted intensity of both laser 
fields. A schematic of the experiment is given in Fig. ^4] 




Figure S4: Cavity optomechanical probing and detection setup. Two orthogonally polarized beams, laser 
1 (dashed purple) and laser 2 (solid red) propagate through the system. Boxes represent polarizing beam 
splitters. Dashed lines indicate optical loss ports. Photodetectors are in blue. Various operators are labeled 
in accordance with the text. 

S.l.l Equations of Motion 

We start from a Hamiltonian H = Hq+H k +Hy describing the mechanical and optical evolution, dissipation 
and interactions [13] : 

Hq = TiuJraC^ c + huj c a[ai + fkj c a2<22 + hGiZ zp (c + c^)a\a\ + hG2Z zp (c + c^)a\a2 (SI) 

where uo m is the mechanical resonance frequency, c is the annihilation operator for the mechanical resonator, 
uj c is the cavity resonance frequency, and ai, a2 are annihilation operators for the two polarizations of the 
cavity mode. Gi, G2 are the optomechanical coupling constants, Z zp = ^jTijlmuJm is the oscillator zero point 
motion, where m is resonator's effective mass. Single photon optomechanical coupling rates are gi = GiZ zp , 
92 = G2Z zp . H K represents the optical input and output coupling, and Hr represents the coupling to the 
mechanical bath. We will eventually specify the results of our calculation to identify the parameters of laser 
1 with the signal beam (5) and laser 2 with the meter beam (m) of the experiment. The Heisenberg Langevin 
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equations of motion are: 




ai{t) = ~ [ai(t),H ] - ^ai(t) + ^/f^La inL1 (t) 

0>2(t) = ~ [a 2 (t),i?o] - ^ a 2(t) + V^ a inL2(^) 

c(t) = -l[ c (t),tf ] + V^/W (S2) 

We have introduced a Langevin noise operator 77(f) to model the thermal mechanical bath, with mechanical 
decay rate IV Additionally, we include input optical field operators for the three coupling channels of each 
cavity mode, input mirror kl, internal loss «i nt , and output mirror kr. We assume the cavity couplings are 
identical for both optical modes. By introducing the definitions for the optical field operators we find: 

ai(t) = (ai + d 1 (t))e~ luJlt , a 2 (t) = (a 2 + d 2 (t)) e -^ t 

flinLi W = feni + £li W + d Xl {t))e l ^\ a inL2 (t) = (a m2 + £ L2 W)e tW2 ' 

flinRi W = 6u Me™ 1 *, a inR2 (t) = ^ 2 (t)e^ 2t 

ainti(t) = eintiWe*" 1 *, a int2 (t) = 6nt 2 (t)e^ 2t 

We allow for coherent field inputs from the left with real valued amplitudes a^ n i = (a^ n i), a^ n2 = (a^ n2 ) and 
frequencies cji, co> 2 , as well as vacuum noise inputs on all ports from the Langevin noise operators £, labeled 
by subscripts indicating the port and laser mode. Classical input intensity noise on laser 1 is modeled with 
the real valued classical noise operator dx±, which is measured in units relative to the shot noise level. d±, 
a 2 represent the complex valued classical amplitude of the intracavity field modes, while the operators di, 
d 2 represent the small classical and quantum fluctuations of the intracavity field, z = Z zp (c + c^) — z is the 
displacement operator, defined so the mechanical coordinate is centered around the optomechanically shifted 
equilibrium location z. Ai = uo c — uj\ + G\z and A 2 = uo c — uo 2 + G 2 z represent the input laser detuning from 
the optomechanically shifted cavity resonance. The equations of motion are linearized by neglecting small 
terms of order d 2 , d x z. To solve the equations of motion we Fourier transform the Eqns. according to 
the following conventions: f(u) = f^e^ f{t)dt, f\uj) = e«" /t (t)dt, {P^rf = f(-u). 



d 1 (uj) = x&{u) (-iGiaiz(u) + Ci (<*>)) 
d 2 (u) = Xc2(^) (-iG 2 a 2 z(uj) + C2M) 
z(cj) __ 1 



2uj m (a^iXciMGM +ai^ix*i(-^)Ci f (^) + ^£2X02(^X2 (w) 
+a 2 ^ 2 x c2 (-a;)C 2 (^)J 



Xm(-^) XmM/ y 

We have introduced the mechanical susceptibility Xm(v) = (To/2 — — cj m )) _1 , and the cavity suscep- 
tibilities for the two modes Xci(^) = (ft/2 — + Ai)) -1 and Xc2(<^) = (ft/2 — + A 2 )) _1 . The op- 
tical noise operators are lumped into Ci(^) = V / ^( < 6i( Co + dxi(u)) + ^/kr£,ri(oj) + \/^int6nti(^) and 
C2M = y/KL£,L2{u) + y/^R^i^) + V^nt^^M- The function A/*(cj) = 

(XmHXm(~ w )) _1 — 2zcj m Z2 P (|ai| 2 GiIIi(a;) + |a 2 | 2 G2lI 2 (cj)) modifies the mechanical susceptibility and 
III (a;) = XciM - Xti(-u), n 2 (cj) = x<a(w) ~ x£ 2 (-<*>)■ 
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The two-sided mechanical displacement spectrum S z (w) = (z(-uj)z(u)) may then be calculated from 
the expectation values of the operator pairs: 

(£li(-^)£Lm) = (6nti(-^)eLi(^)) = (^i(-^)eLM) = i 

(eL 2 (-^)eLM) = (6nt2(-^)eL 2 (^)) = (e^(-^)eLM) = 1 

(^(-a;)^ 1 "^)) = n th + 1, (rf (-u)rj(u)) = n th 
(dxi(-uj)dxi(uj)} — B\ 

Here n t h is the thermal occupation of the mechanical oscillator, and all other expectation values of products 
of Langevin operators are zero. The classical intensity noise is assumed to be locally white and takes a value 
B\ times that of shot noise. Here we assume the Langevin operators as well as the displacement spectrum are 
delta function correlated i.e. (£(— = S(u — uo') with assumed integration over uj' for the calculation 
of all experimentally relevant quantities. 

S?\u) 1 / / n th + 1 n th \ 

+ 4a;^^|ai^iXci(-o;)| 2 +4^^|a 2 ^2Xc2(-^)| 2 

+ 4^^|a m ( X ciM + Xc*i(-^))| 2 ^ij (S3) 

The displacement spectrum consists of four terms. The first term represents the residual thermal motion of 
the optically cooled oscillator. The second term is the displacement due to RPSN from laser 1. Assuming 
laser 2 is responsible for the majority of the optical damping, the third term contains most of the oscillator's 
zero point motion as well as the small RPSN effect from laser 2. The last term is the response to classical 
intensity noise on laser 1. 

Given the classical and quantum fluctuations in the optical spectrum, we want to calculate the mechanical 
response. Let us assume Ai ~ 0, as is the case in the actual experiment. Then A sn = (^^a*a) _1 is the 
output relative intensity spectrum of laser 1 due to shot noise, and A cn = kl\Xci(w) + Xci(~ uo)\ 2 B\/ (a*a) 
is the output relative intensity spectrum of laser 1 due to classical intensity noise. The classical output 
noise reflects filtering by the Lorentzian cavity response, whereas full shot noise appears on the output light. 
However, inside the cavity the shot noise intensity fluctuations are suppressed by the cavity Lorentzian 
[25] . Thus we must treat the perceived level of classical and shot noise differently to correctly infer the 
mechanical response. Using Eq. §we find that Sf /S™ = kk r \x c i(-u)\ 2 ] A sn / 'A cn . Here S s z n and Sf 1 are 
the contribution to the displacement spectrum from the shot and classical noise on laser 1, equal to the 
second and fourth terms of Eq. ^3] respectively. Experimentally this means that even when the measured 
classical intensity noise at an output photodetector is only a few percent of the shot noise level, it may still 
represent a significant amount of radiation pressure drive compared to the shot noise drive. 

S.1.2 Correlation Spectrum 

(2) 

We next turn to the computation of the photocurrent cross correlation spectrum, Sj' (^) = 

((h(— w) — h) — h)) s : with the mean photocurrents I\ = I2 = (l2(t)). Since the pho- 

tocurrents are classical commuting variables, it should be that I\(uS) x 12(00) = h(oo) x I\(uo). To ensure 
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this classical property, we compute the symmetrized expectation value for S/ 12 , defining (f(uj)g(uj)} s = 
(1/2) ((f(uj)g(uj)) + (g(u)f(u))). For comparison with experimental data we compute the one-sided power 

(2) (2) 

spectrum Si 12 (u) = Sj J (lj) + a;). The photocurrents are given by 

h(t) = ^i^i4 etl (t)a det i(t) + Idi(t), h(t) = huj 2 lZ 2 a\ et2 (t)a det2 (t) + J d2 (*) 

where J^i, 7^2 are photodetector dark currents, TZ± = q e ei/TiuJi, 1Z 2 = q e ^2/^2 are the photodetector 
sensitivities (with g e the photoelectron charge) and ei, e 2 include the photodetector quantum efficiencies 
and propagation losses outside of the cavity, adeti and a det2 are the photon annihilation operators at the 
photodetector (see Figure ^4|. The fields at the detector are a combination of the transmitted cavity fields, 
a outi? a out2 and vacuum noise from the optical loss channels. Using the definitions a out i(t) = (a ou ti + 
douti(t))e luJlt , CL ou t2(t) = («out2 + <^out2(^))e* u;2 ^ to distinguish the small fluctuations from the large classical 
amplitude a ou ti = («outi(^)), a out2 = (a out2 (t)), we find: 

a Lti W a deti(£) = eia* utl a out i + 6i(a* utl d ou ti 

a det2(*) a det2(i) = e2a* ut2 a ou tl + e2(aout2^outl (*) + a out2 dl utl (t)) + Ve2(l-e2)(^ut2^n2W + «out 2 ei 2 W) 
where £ n i, £ n2 are Langevin vacuum noise operators. Substituting in the above relations, we find: 

ff/iaM = / ( ^outWoutl(-^) + a utl<Ll(-^A f ^out2^out2(^) + goutgjLtg M 
^1^2 \ \ ^outl^outl y y «out2^out2 

where we have employed the fact that J^i, ^d2, £ni, £n2 are all uncorrelated with each other. Thus terms 

f2) — 

proport iona 1 to expectation values of products of these operators are zero, and Sj i2 (oj) / I\I 2 becomes inde- 
pendent of ei, e 2 . We note, however, the power spectra of the individual photocurrents will depend on J^i, 
^d2, ei, 62, so the ratio Sj 12 / Si 1 Si 2 does improve with increasing detection efficiency and a lower photodetec- 
tor noise floor. We can relate the intracavity photon operators to the output operators using the boundary 
conditions a ou ti = y/KRa,i, «out2 = y^Ra 2 , and d ou ti = V^R^i ~ 6*1, d out2 = ^/^Rd 2 - £r 2 . Applying the 
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solution to the equations of motion from above, we find: 



hh 



l 



I^rGi 

KrGi 

^rGi 



KR\Q>ir\Q>2 

^/K^G 1 \a 1 
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a 2 
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a 2 
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2 

2 zIIi 

2 zIIi 

2 zni 
2 zni 

2 xii 1 



ni(-^)v^G 2 |a 2 | 2 n 2 (cj) (z(-u)z{u)) a 

Uj)a 2 ^KLKRXc2(u) (z(-Uj)£ L2 (uj)} s 

uj)alyjK J - int K jR Xc2{u) (z(-uj)^ nt2 (u)) s 
uj)a 2 (K R Xc2(u) - 1) (z(-u)€ R2 (uj)) a 
uj)a 2 ^K L K R x* c2 (-uj) (z(-uj)£} L2 (uj)^ 

a;)a 2 (^Xc 2 (-^) - 1) (^(-^)£LM) s 

VklKrXcI (-^) ((Cli(-^)^H) s + (^i(-cj)z(cj))J 

u)o>i(krXci(-u) - 1) (^i(-cj)^(o;)) s 
(j)aiV^L^RXciH {(^li(- uj ) z ( u; )) s + (^i(-^)^H)s) 



2 *n 
2 *n : 
2 *n 



2 *n 



^n 2 



2 zn 2 



(j)ai(^x*i(^) - 1) (^i(- CJ ) 2: ( CJ )) £ 



(S4) 



The necessary operator expectation values are: 



(z(-uj)^ j2 (uj)) s = 



g 2 Z zp a 2 x* c2 (u), for j = L, hit, R 



(^ jl (-uj)z(uj)) s = 



(dx 1 (-uj)z(uj)} 6 



#iZ zp aix*i(-^), for j = L,int,i? 
^l^zp (a*Xci(u) + aix*2(-^)) 



To understand which terms are relevant in the experiment for RPSN, we note that IIi(a;) approaches zero 
as Ai approaches zero. The first seven terms of Eq. ^4] then vanish under these conditions. This leaves 
terms that are proportional to correlations between z, the resonator position, and vacuum noise operators 
£ri which represent the shot noise optical driving force on laser 1. In the limit where Ai = and in the 

absence of classical intensity noise, the cross correlation takes a simple form -7^7^ = ^Wi92^m U2 tt^j ~ • 

For classical intensity noise an additional term is added: 2zgig 2 (jj m KL U2 ^^ Xcl ^^ cl ^~^^ B\. Note that 
in our case where T m <C ft, the classical and shot noise terms show similar functional form near uj = uj m . 
Importantly, the shot noise driven term includes an extra phase shift of Arg(x* 1 (cj m )) ~ arctan(2d; m /ft). 
This additional phase shift allows one to experimentally distinguish the quantum versus classical origin of a 
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radiation pressure drive. As shown in Figure £[5)3, the phase offset in the cross correlation varies continuously 
as the classical noise level is increased relative to the shot noise level. 

When Ai is nonzero, the first term of Eq. ^produces a correlation induced by the mechanical motion 
of the resonator imprinted onto both photocurrents. Depending on the sign of Ai the correlation from 
thermal motion may add either constructively or destructively with the RPSN correlation, leading to a laser 
frequency dependent lineshape. Example expected cross correlation lineshapes are presented in Fig. ^5}^. 
Note that this effect is most pronounced near the mechanical resonance peak. Over a wide range of Ai the 
Lorentzian wings of the mechanical resonance are insensitive to the thermal motion and give an accurate 
representation of the RPSN correlation. 




1.54 1.55 1.56 1.54 1.55 1.56 

Frequency (MHz) Frequency (MHz) 



Figure S5: Cross correlation spectra. The magnitude squared (A) and the phase (B) of Si 12 (described as 
Si SM in the main text) are plotted for the parameters of device B. In (A) Ai is varied. The tallest curve 
corresponds to Ai/27r = 10 kHz, and the lowest curve to -10 kHz. The black curve corresponds to the 0.28 
kHz, the value of the correlation measurement in Fig 3A. In (A) the level of classical intensity noise relative 
to shot noise on laser 1 is varied. The bottom curve is pure quantum noise and the top curve is pure classical 
noise. The black curve corresponds to a measured noise at the photodetector consisting of 91% shot noise 
(75% of radiation pressured drive from shot noise) in the correlation measurement shown in Fig. 3A. 



S.1.3 Optical Bistability 

One limit to the amount of optical power circulating in an optomechanical system is the onset of optomechan- 
ical bistability. Here, the mean radiation pressure on the mechanical element causes a static displacement, 
which shifts the cavity resonance frequency by ~ G\z. When this shift becomes comparable to the cavity 
linewidth, there exist two stable values of the circulating power and z for a given input optical power, over 
some range of laser frequencies where the laser is nearly resonant with the cavity. In this situation the 
system may be driven by small noise sources from one stable equilibrium to the the other in an uncontrolled 
fashion. To some extent, such processes can be avoided by using active feedback, as is employed in the 
actual experiment, increasing the duration of stable operation. To estimate the threshold for optomechan- 
ical bistability, we need to consider the net displacement of all the mechanical modes of our system with 
appreciable optomechanical coupling, not just the mode of interest for RPSN measurements. The critical 
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value of the intracavity photon occupation, Nf , above which bistable behavior is observed may be estimated 
for our membrane system by: 

m,n (m,n) 

Here the sum is over mechanical eigenmode indicies (m, n), and tj( m>n ), G( m?n ) are the mechanical resonance 
frequency and optomechanical coupling constant for the mode (m,n). We estimate G( m?n ) by measuring 
the location of the optical mode spot on the membrane, and then computing overlap integrals between the 
expected optical and mechanical mode functions for each mechanical mode. For our measured parameters 
we obtain N{ ~ 3.5 x 10 8 which is comparable to the highest value of Ns employed in the actual experiment. 



S.2 Experimental Methods 

Here we describe in detail the calibration and operation of our optomechanical system. More details about 
the construction of the system can be found in Ref. [24] . 



S.2.1 Experimental setup 

The optical cavity consist of two mirrors, one flat, one with a 5 cm radius of curvature, both with approx- 
imately 1 x 10 -4 fractional intensity transmission at the operating wavelength of 1064 nm and only a few 
xlO -6 scattering and absorption losses. The mirrors are held at a separation L = 5.1 mm by an invar spacer. 
With this geometry, the cavity is expected to support a Gaussian profile standing wave mode with a 1/e 2 
intensity radius of 72 /am. The flat mirror is attached to the invar spacer through a multilayer pizeoactuator 
that is used to precisely tune the overall optical path length of the cavity. The empty cavity finesse is 
measured to be 31,000 via optical ringdown spectroscopy. 

A stoichiometric high-tensile-stress SiaN4 membrane supplied by Norcada Inc. is placed inside the cavity 
as the mechanical element. The membrane is square in shape, 0.5 mm on a side and is suspended on a silicon 
frame with dimensions 5 mm x 5 mm x 0.5 mm. The membrane film is typically measured via ellipsometry 
to have a refractive index of 2.0 and film thickness of 40 nm, yielding a reflectivity of 10 % for 1064 nm light. 
The mechanical modes that most strongly couple to the light are transverse "drumhead" modes, which in 
the high-tension limit, show sinousoidal displacement profiles, with motion out of the plane of membrane. 
We label the modes with two indices (i,j) denoting the number of antinodes of oscillation along each of the 
transverse direction. For most of the experiments performed we focus on the (2,2) mode which has four 
antinodes, one in each corner of the membrane. This mode oscillates at a frequency uj m /2iT ~ 1.55 MHz, 
which is twice the frequency of the fundamental (1,1) mode. The intrinsic mechanical linewidth, which varies 
with temperature and mounting technique, is measured in-situ, at a wavelength where the cavity finesse is 
low to avoid any optomechanical effects, via mechanical ringdown spectroscopy, giving Tq/2k = 0.116 Hz for 
device A and Tq/2tt = 0.47 Hz for device B. 

The membrane is positioned in the cavity, so that the optical mode spot is approximately aligned with 
one of the antinodes of the (2,2) membrane mode to attain the largest optomechanical coupling. The 
membrane is located about 0.9 mm from the flat mirror, and its position can be finely tuned along the 
optical standing wave with another multilayer piezoactuator. The effective cavity input and output coupling 
and loss as well as the optomechanical coupling all vary with the location of the membrane along the optical 
axis. In conjunction with a simple matrix model [23], we use optical ringdown measurements of the cavity 
linewidth when the membrane is placed at the operating point to extract values for the cavity parameters. 
We determine n the total cavity linewidth, and k,r, the input and output coupling of the cavity, and 
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ftint, the internal cavity loss coupling. fti nt is due mainly to membrane absorption and scattering as well as 
clipping of the optical mode by the membrane frame. This analysis yields ft/27r=1.17 MHz, ftz, = 0.33ft, 
kr = 0.23ft, and fti nt = 0.44ft for device A and ft/27r=0.89 MHz, kl = 0.32ft, kr = 0.59ft, and fti nt = 0.09ft 
for device B. 

A detailed diagram of the laser setup employed in the experiment is shown in Figure A diode- 
pumped, monolithic non-planar ring oscillator type Nd:YAG laser from Innolight GmbH drives the entire 
experiment. The laser is first spectrally filtered by passing through a 40 kHz Fabry- Perot optical cavity. 
Next, the light is doubled passed through an acousto-optical modulator (AOM) used for fast control of 
the laser frequency. The laser is then split into two components that serve as the signal and meter beams. 
Each component is passed through another AOM, which provides relative frequency control and independent 
intensity stabilization of the two beams. The signal beam passes through an in- fiber electro-optical modulator 
(EOM) that adds frequency sidebands at 18 MHz for Pound-Drever-Hall frequency stabilization of the laser- 
cavity detuning. At low frequencies the overall optical path length of the cavity is servoed to maintain the 
cavity resonance with the signal beam. At higher frequencies, up to 100 kHz, the laser frequency is servoed 
via the common AOM. With this servo system we are able to reduce the cavity-laser frequency fluctuations 
to the few kilohertz level. The servo output is aggressively filtered to ensure there is no response near 
uj m . The location of the membrane is passively stable at the few tens of nanometers level during cryogenic 
operation. To combat drift in the servo lock point, residual amplitude modulation from polarization drift 
in the EOM is actively canceled. To accomplish this we apply a DC voltage to the EOM crystal to null the 
amplitude modulation measured by a photodetector sampling the beam after the EOM. With this system, 
we typically see frequency offset drifts of less than 1 kHz over the tens of minutes times scale relevant to 
data taking. The signal and meter beams are combined on a polarizing beam splitter, and their polarizations 
are rotated to match the polarization eigen-axes of the slightly birefringent cavity (~ 400 kHz birefringence 
splitting). After the cavity, the two beams are split by another polarizing beam splitter and directed on 
to individual photodetectors. We typically see less than 10 -3 cross coupling between the two beams. The 
photocurrents are simultaneously digitized with a two-channel, 16-bit digital oscilloscope card and a discrete 
Fourier transformation is applied to the digitized signals to generate power spectra and cross correlations. 

The optomechanical system is mounted onto the cold finger of a low- vibration 4 He flow cryostat from 
Advanced Research Systems Inc., which attains a base temperature of 4.9 K. The cryostat has windows that 
allow direct free-space optical access to the cold sample region. We include an additional copper radiation 
shield with small-aperture, fused-silica windows around our system, heat sunk to the cold finger to ensure 
our device thermalizes to the cryostat base temperature. In cooling from room temperature to cryogenic 
temperatures the membrane retains sub-milliradian alignment with the cavity optical axis. We are not 
inhibited by vibrations from cryogen flow in our system in part due to the low-vibration cryostat design 
made to eliminate boil off near the cold finger, and in part due to the near-monolithic design of our cavity 
[24] and high active feedback bandwidth. 

S.2.2 Calibration 

The optomechanical coupling of system is calibrated in two ways [24]. Because the signal and meter beams 
have the same spatial profile, we assume they share a common single-photon optomechanical coupling rate 
g = gs = gjw and coupling constant G = gZ zp . At low signal laser power, the effect of the meter beam on 
the resonator can be predicted by calculations similar to those in the theory section, as is well documented 
in e.g. Ref. [25 . In the limit where To, <C Tm, the mechanical damping is given by T m = Tm = 
g 2 NMK (\Xc2(^m)\ 2 — \Xc2(— ^m)| 2 ) • Using this relation we can extract g from r m |Ar s =o and a measurement 
of Nm = Im / (Qe£2KR) and the cavity parameters. With this method we estimate g/2ir = 15.7 Hz for device 
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Figure S6: Detailed optical setup. The signal beam (dashed purple) and meter beam (solid red) are derived 
from a single passively filter 1064 nm source. Acousto-optical modulators (AOM) are used to shift the laser 
frequency. An electro-optical modulator (EOM) is used to apply frequency sidebands for a Pound-Drever- 
Hall (PDH) frequency lock. Light is detected at various points using photodetectors (PD). Dashed black 
lines represent beam splitters, boxes represent polarizing beam splitters. A/4 and A/2 are quarter- wave and 
half-wave retarders respectively. 



A and g/2it = 16.4 Hz for device B. Alternatively, we can use the thermal motion of membrane as a known 
displacement to calibrate g. The effective temperature of the membrane mode is expected to be given by 
T e ff = TbathPo/r m . We measure Tb a th from the silicon diode thermometer on the cryostat to be about 4.9 
K. Then g 2 = {S lM \ Ns=0 {u m )/P M ) Y 2 m hw m / (8|n 2 (cj m )| 2 ^T bath r ) , yielding g/2n = 16.9 Hz for device 
A and g/2ir = 15.8 Hz for device B. For the cross correlation measurements on device B, we measure a 
slightly smaller mechanical damping rate for the same laser parameters as compared the data in Fig. ID. 
Thus we estimate a slightly smaller value of g/2it = 14.8 Hz for this data. This difference can be accounted 
for by a small drift in the membrane position along the optical standing wave gradient between the two 
measurements. 

We estimate the detuning of the signal beam in a variety of ways. For the data of Fig. ID, we look 
for linear trends in both T m and the mechanical resonance frequency (the so called "optical spring" effect) 
as a function of Ns- Using these trends, we estimate As/2ir = 1.5 kHz for device A and As/2tt = 2 kHz 
for device B. Alternatively, using the data of Fig. 3 A, we fit the cross correlation to the functional form 
derived above, with A^ as a free parameter. This fit yields As/2tt = 300 ± 100 Hz. As mentioned above, the 
small drift in A^ over the time between the measurements accounts for the difference in the estimates for 
device B. Below, in the analysis section, we show how to account for this coherent cooling or heating from 
the signal laser in order to accurately assess the contribution of RPSN from the mechanical displacement 
spectrum data. Additionally, as discussed in the theory section, with nonzero A^, the cross correlation 
becomes distorted and contaminated by thermally driven mechanical motion. 
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The measured classical intensity noise on the signal beam after the cavity is always a small fraction of 
shot noise. We assess its value and impact on the experiment in several ways. Direct measurements of 
the total photocurrent power spectrum contain both shot noise and classical noise contributions, as well 
noise from photodetector dark current. We calibrate the shot noise plus detector noise by illuminating the 
photodetector with a shot noise limited incandescent light source producing the same average photocurrent 
as the laser light used in the experiment. We subtract the resulting photocurrent spectrum from Si s (uj) 
measured during the experiment, and expect the residual to be the classical noise. The recorded shot noise 
spectral density is within a few percent of the predicted value, 2q e Is: based on a careful measurement of 
the photodetector transimpedence and digitization electronics frequency-dependent gain. However, as the 
classical noise is typically less than 10% of the shot noise level, this method is subject to significant error. 
From these measurements we extract a classical noise power spectral density of less than -157 dBc/Hz. We 
believe the majority of this noise arises from intensity noise imprinted on the laser from the electronics that 
drive the second AOM on the signal beam path as shown in Figure Because the meter beam passes 
through an independent AOM with independent electronics, we do not expect any classical noise on the two 
lasers to be correlated. Thermally occupied mechanical modes of the mirrors and cavity support structure 
lead to noise in the cavity resonance frequency in the megahertz range. However, because the signal beam 
is operated close to the cavity resonance, we do not expect a large effect from frequency to amplitude noise 
conversion by the cavity. We experimentally confirm that we are insensitive to such effects, because our 
signal beam noise floor is, to good approximation, independent of As. 

Classical laser intensity noise can also be detected and differentiated from shot noise by its effects on 
the membrane mechanics. As pointed out in the theory section above, the phase of the cross correlation 
spectrum depends on classical versus quantum noise drive. By fitting the phase of Si SM (cj), we estimate that 
75% of the correlation signal is due to RPSN, with most of the remainder accounted for by classical radiation 
pressure noise drive. This value is within the range of the estimate made by direct photodetection of the 
classical noise. Moreover, by taking into account the measured 25 degree phase shift in the photodetection 
electronics at frequencies near the mechanical resonance, we find excellent agreement between the absolute 
phase of the black theoretical curve of Fig. ^5^B) and the red measured curve of Fig. 3C. 

Another measure of the classical noise level can be extracted from the RPSN data of Fig. ID (reproduced 
in Fig. ^7|). Here, we expect the increase in peak spectral density due to RPSN to scale linearly with the 
intracavity photon number. However, for a constant classical relative intensity noise, the peak spectral 
density should scale quadratically with the intracavity photon number. A three term polynomial fit to the 
data yields a constant term that is related to the thermal motion, a linear term related to RPSN, and a 
quadratic term related to classical radiation pressure. This method indicates that about 75% of the increase 
in peak spectral density is due to RPSN in device B for the highest N$ value employed. 

At large circulating optical power, it is possible that absorbed laser light will heat the membrane and 
increase Tb a th- This effect, in the simplest case, would produce a linear increase in the peak displacement 
spectral density as function of optical power, the same scaling as RPSN. To control for laser absorption 
heating, we look at the (4,4) membrane mechanical mode. For the (4,4) mode the optomechanical coupling 
is weaker, the environmental coupling is stronger and the oscillation frequency higher. Thus we expect only 
few percent increase in the peak displacement spectral density due to RPSN at the highest drive power. Any 
significant increase in the amplitude of motion of the (4,4) mode, we may then attribute to spurious heating 
effects. The data of Figure ^7|show the response of the (2,2) and (4,4) mechanical modes for both devices. 
There is no clear linear trend in either of the (4,4) mode data sets. This allows us to put an upper bound of 
10 % for the increase in bath temperature at the largest N$ employed for both devices, if we assume all of 
the mechanical modes are coupled to a common thermal bath. 
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Figure S7: Response of multiple mechanical modes. The peak spectral density is plotted for: part (A) device 
A (4,4) mode (blue diamonds), part (B) device A (2,2) mode (yellow squares), part (C) device B (4,4) mode 
(red diamonds), and part (D) device B (2,2) mode (green squares). The (2,2) modes of both devices show a 
strong linear trend in response to RPSN, while the (4,4) modes show a relatively weak response. 



S.2.3 Data and Analysis 

For each point of the device A data in Fig. ID, we average the power spectra of 1000 records of the 
photocurrents generated by signal and meter beams. Each record is 20 ms in length and is digitally sampled 
at 5 x 10 7 samples per second for each laser power setting. We compute Si m (uj)/I m and convert this into an 

effective displacement spectrum using the relation: S z (uo) = Sl f 2 ^ — — 1 . Over a region encompassing 
several mechanical linewidths around the mechanical resonance frequency, we fit each curve to a Lorentzian 
profile to extract the mechanical linewidth and peak spectral density. For each value of Ns we record the 
response of the system for several values of the meter beam power, corresponding to mechanical damping rates 
ranging from 0.8 to 4 kHz. We perform a linear fit to T m versus Nm and fit S z (uj m ) inversely proportional 
to Nm, to extract an estimate of S z (uj m ) at r m /27r = 3 kHz. 

For the device B data, in Fig. ID, the signal to noise is much greater, so we only need to average a few 
hundred records for each laser power. We use only a single meter beam power for each value of N$. Again 
we fit the averaged spectrum to a Lorentzian to extract T m and S z (u m ). Because A s is a few kilohertz for 
this data, it provides an additional optomechanical damping and T m shows a linear trend for increasing Ns, 
with its value decreasing by about 20 % at the highest signal beam power. By extrapolating to Ns = 0, we 
find T m approaches 2tt x 1.43 kHz in the absence of the signal beam. We then apply a small correction to 
the S z (u) data to remove the effect of the optical damping from signal beam. 

The dashed theory bounds of Fig. IB are generated using the calibration methods discussed above. 
The boundaries are generated by finding the extrema of the theoretical estimate of S z (uj m ) using the two 
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estimates for g and the range of estimates for the classical intensity noise contribution. The horizontal error 
bars on the data represent the systematic uncertainty in the conversion of mean photocurrent to N$. The 
vertical error bars represent statistical error in the measurement. 

For the cross correlation data we take 1000 records of 20 ms in length digitally sampled at 2 x 10 7 
samples per second. We simultaneously record both photocurrents using a two channel digital oscilloscope 
card. We use vector averaging to compute the complex cross correlation spectrum, and from the same 
data also compute the scalar average of the power spectra of the two individual photocurrents. We believe 
that 1000 averages is sufficient to converge the cross correlation spectrum in region of several mechanical 
linewidths in spectral width. As evidence, the blue data in Fig. 3B has converged to the same lineshape as 
the data in Fig. 3A, despite the added mechanical noise. The boundary of the gray theory band of Fig. 3A 
is again generated as the extrema of the theoretical estimate using estimates for g and the classical intensity 
noise level. 
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